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We study charge transport of a topological superconductor connected to different electromagnetic 
environments using a low-energy description where only the Majorana bound states in the super¬ 
conductor are included. Extending earlier findings who found a crossover between perfect Andreev 
reflection with conductance 2e^//i to a regime with blocked transport when the resistance of the 
environment is larger than 2e^/fe, we consider Majorana bound states coupled to metallic dots. In 
particular, we study two topological superconducting leads connected by a metallic quantum dot in 
both the weak tunneling and strong tunneling regimes. For weak tunneling, we project onto the most 
relevant charge states. For strong tunneling, we start from the Andreev fixed point and integrate out 
charge fluctuations which gives an effective low-energy model for the non-perturbative gate-voltage 
modulated cotunneling current. In both regimes and in contrast to cotunneling with normal leads, 
the conductance is temperature independent because of the resonant Andreev reflections, which are 
included to all orders. 


I. INTRODUCTION 

There is currently a large attention towards topological 
superconductors, which are expected to host zero-energy 
states with interesting properties. These states, known as 
Majorana bound states (MBSs), constitute half fermions 
in the sense that two MBSs are needed to define a sin¬ 
gle fermionic level, where the degree of freedom of this 
fermion corresponds to the parity of the fermion number. 
Thus, if the MBSs are well separated in space the parity 
degree of freedom is a topologically protected quantity, 
which can potentially be used for quantum computation 

m- 

Promising candidate systems are hybrid materials 
where s-wave superconductors are used to proximitizes 
other materials with strong spin-orbit coupling US]. 
Systems like this have already been studied in a num¬ 
ber of experiments ISHS]. 

An important diagnostic tools to verify the presence of 
MBSs is tunnel spectroscopy for which the zero-energy 
bound states are predicted to give resonant Andreev 
reflection and hence a low-temperature conductance of 
2ev/i mn]. With interactions in the normal leads the 
situation changes because the Andreev reflection may be 
suppressed by repulsive interactions in the metal. This 
was studied by Fidkowski et al. who showed that 

when the normal metal is a Luttinger liquid, the resonant 
Andreev reflection is replaced by a non-conducting fixed 
point for interactions stronger than a critical value. In 
particular, the prediction is that for a chiral Luttinger 
liquid the cross-over happens for the Luttinger parame¬ 
ter AT = i, with the insulating phase occurring for strong 
interactions, AT < |. More details of the voltage and tem¬ 
perature dependence were studied later m- A related 
study of a point contact between a Luttinger liquid and 
a chiral Majorana mode was done in Ref. [HI 

Another interesting situation is when the Majorana 
bound state is tunnel coupled to an interacting quan- 
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FIG. 1. Illustration of the three types of junctions consid¬ 
ered here. In (a) there is a junction between a topological 
superconductor (SC) and metallic lead with its electromag¬ 
netic environment modeled by an impedance Zenv(u)). In (b) 
the environment is replaced by a metallic island connected to 
the drain by a quantum point contact. In (c) is a junction 
between two topological superconductors, where the junction 
consists of a metallic Coulomb blockaded island. 


turn dot [TMT!?] . For such a setup the Coulomb block¬ 
ade results in sharper resonance structures in the weak- 
coupling limit |15j . This was shown to hold true also for 
stronger coupling by Cheng et al. [T5], who showed that 
the perfect Andreev reflection dominates over the Kondo 
effect in the strong-coupling hxed point. 

In this paper, we analyze the transport proper¬ 
ties through MBSs via different types of environments, 
namely (a) a general linear electromagnetic environment, 
(b) a metallic quantum dot with a quantum point contact 
drain, and (c) a metallic dot with a second topological 
superconductor as a drain. See Fig. for illustrations. 
The situation in (a) was shown by Liu |20] to exhibit a 
crossover from conducting to non-conducting when the 
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environment impedance is larger than the resistance of 
the resonant Andreev conductance > /i/2e^, and 

the same results follows in the situation (b) when the 
quantum point contact has less than two open channels. 
The nature of this transition is completely analogous to 
the Luttinger-liquid physics studied by Fidkowski et al. 

m- 


The above transitions occur when the isolated MBS- 
metal junction resistance and the environment resistance 
are equal. This motivates the study of situation (c) in 
Fig. 1 where two MBS-metal junctions are connected via 
a Coulomb-blocked metallic island. Because both junc¬ 
tions have resistance /i/ 2 e^, the combined resistance is 
locked to the transition value and the conductance is 
temperature independent, similar to a non-interacting 
scattering problem. In fact, after integration over charge 
fluctuations the resulting low-energy model can be solved 
by mapping to a simple fermionic scattering problem by 
refermionization. 

The paper is organized as follows. In Sec. [llj we in¬ 
troduce the key elements of the models, namely a single 
MBS coupled to a fermion bath and the bosonic electro¬ 
magnetic environment, and show how these elements are 
bosonized. In Sec. |III[ the effective low-energy actions 
is derived and used to identify the crossover impedance. 
The tunneling limit is treated in detail in Sec. |IIIB| The 
case of a quantum point contact is included as special 
case in Sec. Ill C| Next, the Coulomb blockade of a metal¬ 
lic island coupled to two MBSs is analyzed in Sec. [Tv| 
using the dual bosonized representation. Finally, conclu¬ 
sions are given in Sec. |V] 


II. BOSONIZATION OF THE MBS-NORMAL 
METAL-ENVIRONMENT CONTACT 

We start with the Hamiltonian of a MBS tunnel cou¬ 
pled to a normal-metal lead, with the inclusion of an dis¬ 
placement of the environment charge which accompanies 
the tunneling event 

H = Hgi + I?env + .ffx.env; (1) 

where 

-^ei (2) 

k 

i?T.env =t ^ 7, (3) 

k 


Because the MBS couples to a particular spin projection 
of the lead m , we only need to include the corresponding 
spin direction and we can effectively use a spinless model. 

Moreover, we take the tunneling matrix element t to be 
independent of k (s-wave scattering) and, furthermore, 
assume that the density of states is constant: p(^) = 
<5 (C ~ Cfe) = Po- These simplifications allows us to 
parameterize the dispersion relation as = Vp (k — kp) 
and map i?ei to a ID chiral electron gas 


/ OO 

dx ip^{x){-idx - kF)'ilj{x), (5) 

-OO 


with 


Cfe = 4 ^ [dxe ( 6 ) 


v^, 


and C =2ttpqVp playing the role of a normalization 
length. In this language, the tunneling Hamiltonian be¬ 
comes 


FlT = A(V'(0)e*“-V'^(0)e-“)7, (7) 

where A = -\/27rpo%?. 

This ID Hamiltonian can now be bosonized in the stan¬ 
dard way, for example by folding into a half-infinite wire, 
such that tp{x > 0 ) = iPr{x) and ipi^x < 0 ) = ^^^{x). 
Following Fidkowski et al., we integrate out all x ^ 0 
degrees of freedom, which then gives the action 

S = So[e]+ST[Q,a] + S,^M, 

5'o[0] = H |wn||0(*w„)|2, 

iu)n 

rP 

5't[ 0] = 2A' / dr Ua, cos (0(r) -|- a(T)), 

Jo 

where 0 (r) relates to current at x = 0 in the half¬ 
wire language (Fidkowski et al. [12]) or to 0 = 
TT J_^dx sign(x)p(x) in the ID chiral language. Fur¬ 
thermore, the Pauli matrix cr^ refers to the two-level sys¬ 
tem defined in a complex-fermion basis given by 7 and 
another MBS in the topological SC (for example 72 in 
the Fig. 1(a)). Also a short scale cut-off parameter a 
has been included, so that = Xj^J^'Ka. Similarly, the 
current operator in its bosonized form is given by 


( 8 ) 

(9) 

( 10 ) 


and where k are the quantum numbers for the normal- 
metal eigenstates that couple to the MBS, 7 , and Hem 
is the Hamiltonian of the linear environment. The oper¬ 
ator a is the canonically conjugate operator to the en¬ 
vironment charge Q, i.e. [a, Q] = i. Therefore 
creates/annihilates a charge in the environment respec¬ 
tively. The current operator is given by 

/ = *et^(c,e- + cIe-“)7. (4) 

k 


/(r) = 2A'(Ta; sin (0(r)-I- q;(t)) . (11) 


III. MBS COUPLED TO AN 
ELECTROMAGNETIC ENVIRONMENT 

In this section, we study a MBS coupled to two differ¬ 
ent environments: firstly a general linear impedance and 
then secondly a quantum point contact (QPC). 
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A. MBS coupled to a linear environment 


The environment action describes the dynamics of the 
field a and thus after the integration over all other envi¬ 
ronment degrees of freedom, it can be expressed as 

S'envM = ^^G-^{iuJn)\aiiu}n)\^, (12) 

itOn 


where a’s Green’s function is dehned as Gaa{T) = 
— (TT-a(T)Q;(0)). The a correlation now follows from 
classical circuit theory, analogous to the environmental- 
Coulomb-blockade theory [21}{23]. In frequency domain 
we have 

^ ^ e'^\Zt{iuJn)\'^ 

Gaa(*W„) = ^- ^Gvv{^^n) = - -^-•G//(^W„), 

(13) 

where Gyv and G// are the voltage-voltage and current- 
current correlation function, respectively, which by the 
Kubo formula are related to the total impedance of the 
network as G// = jcUnlZ^^. We thus get an expression 
for Gaa in terms of the total environment impedance 

e^Z^(iuJn) 


Gcta{}^n) — 


\UJn 


(14) 


The network consists of the junction capacitor and the 
environment impedance in parallel (see Fig. 1(a)) and 
therefore we obtain 

GaaiiUn) = -7--TW , (15) 


(K|G + (Z,*„J-y 


where Z^^^ is the impedance of the network connecting 
the junction to the voltage source. 

Now, from the action (10 1 we see that the field 0 — a 


does not couple to the junction and we can integrate it 
out. This leaves us with an effective action for 0+ = 
0 -I- a where the quadratic part becomes 


5o[0+] = ^^K|G;i(*a;„)|0+(zc.„)| = 


27r/3 


(16) 


with 


G+{iuJn) — 1-1- Gaaii^n), 

TT 


(17) 


and the non-quadratic part is given by Eq. (10 1 . 

Let us now study the case where the environment 
impedance can be approximated by a constant real 
impedance. For small frequencies the Green’s function 
G_|_(ia;„) reduces to 


G+(ia;„) « 1-I- 7 , (18) 

where we have defined 7 = 2e^Zenv/^- The environment 
model is thus identical to the situation of a MBS coupled 
to a Luttinger liquid with the replacement K (l-l-q)”^. 
Thus, when the environment impedance exceeds hl2e^, 
corresponding to 7 > 1 and K < 1/2, there is a transition 
from resonant Andreev reflection to insulator behavior at 
low temperature and low voltage. 


B. Line shape of Majorana environmental Coulomb 
blockade close to the non-conducting fixed point 

It is also interesting to discuss the finite-volt age line- 
shape of the peak in differential conductance, because 
this is most often the experimental signature used to 
conclude on the presence of MBSs in the superconduc¬ 
tor. The tunneling current is calculated to second order 
in t by standard perturbation theory: 

I=-iej dt({/(0),7LT.env(t)}). (19) 

J —00 

After some algebra, we then obtain 

/*oo 

I=\] du:[f{uj-eV)-f{uj + eV)]P{u:), (20) 

where f{uj) = I/(exp(/3tLi) -|- 1) is the Fermi-Dirac dis¬ 
tribution function where we sat the reference chemical 
potential of the superconductor to zero, F = 27re^t^p 
(with p being the density of states in the normal-metal 
contact) and where 


P{u:) = 






( 21 ) 


is the usual “P(i?)-function” [^, which obeys 

P{u)) — P(—tLi)e^“ and describes the response of 

the environment. Without coupling to the environ¬ 
ment, we have P(w) = 2'k5{oj) and hence dl/dV = 
47 rF( 2 e^/h)(—/'(eE)). With an environment the differ¬ 
ential conductance is 


dl 

dE 





which at zero temperature reduces to 


eV)), 

( 22 ) 


dJ 

dE 


r=o 


2 e2 

h 


F(P(eE) -f P(-eE)). 


(23) 


Interestingly, the normal-metal-MBS junction thus mea¬ 
sures the P(w), in contrast to the usual metallic junction 
where the differential conductance is given by the inte¬ 
gral of P(w) [H]. At zero temperature and low energies, 
the environment function goes as P{uj) ex w'*', which con¬ 
firms the conclusion from the previous section that for 
7 > 1 , the zero-bias conductance goes to zero at low 
temperatures. 


C. Metallic dot with a quantum point contact 
drain 

Next, we consider a special kind of environment, 
namely a metallic dot with connections to a MBS and 
quantum point contact with N open channels to the drain 
electrode, see Fig. 1(b). The environment action includes 
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the charging energy of the dot, the coupling to the phase 
field a, and the bosonized open channels of the quantum 
point contact. It reads 

fP ( e 2 / 

5'env = J dr I iaQ + — I Q + - 

+ (24) 

i,n 

where Q is the conjugate field to a, counting the num¬ 
ber of charges passing through the Majorana-bound-state 
junction, and 4>i is the boson field describing the charge 
that passed from lead to dot through mode i in the QPC 
[211 US], and Ng is controlled by a gate voltage. After 
integrating out all (piS as well as the field Q (and remove 
Ng by a gauge transformation) the environment action 
describing the field a becomes identical to the electro¬ 
magnetic model above with the environment impedance 
given by Zg^v = h/Ne^ or in terms of the parameter 7 
defined above, we get an effective 7 for the QPC setup 
given by 

7 qpc = 2/iV. (25) 

This shows that the there is a crossover between full nor¬ 
mal reflection and full Andreev reflection at N = 2. 

IV. TWO MBS COUPLED VIA A METALLIC 
DOT 



Finally, we consider the device setup in Fig. 1(c) con¬ 
sisting of a metallic Coulomb-blocked island coupled to 
two MBSs. It is an inverse version of the system studied 
by Fu [21] and later by Hutzen et al. [ 21 HH], namely 
a topological superconductor Coulomb island with the 
two MBSs coupled to normal leads. The setup allows 
for a non-perturbative solution of the conductance and 
therefore gives an interesting way to investigate the com¬ 
bination of MBSs and Coulomb interactions. The Hamil¬ 
tonian of the metallic dot coupled to two MBS is 


H= Y Aei,,±AT,.±Ac, 

r^L,R 

(26) 

Ael,r = ^ 
k 

(27) 

HT,r = tr ^ ( {flk,r ~ 
k 

(28) 

He = Ac(Al ± Ar - Ng)\ 

(29) 


where r = L, R labels the two ends of the quantum dot, 
which have independent quantum numbers labeled by fc, 
and the MBSs are denoted by 7l,r- The current opera¬ 
tors at the two contacts are 

It = ietr -f cl^^ Jr- (30) 

k 


For weak tunneling and away from Ng being a half¬ 
integer, the dot has a well defined charge and current is 
therefore blocked as in the usual Coulomb blockade. In 
this situation, current is carried by cotunneling, which 
is via virtual states with one additional charge. Second- 
order perturbation theory then gives the following result 
for the conductance at T ^ Eq 

Gcotu„ = ^rLrR(^^ + —) , (31) 

where F^ = 27ruF|trP and (1 T 2 A"g), with 

Ng e [—1/2,1/2]. We note that the cotunneling current 
is independent of temperature at low temperatures. The 
cotunneling conductance diverges near the charge degen¬ 
eracy points Ng = ±1/2. To understand the behavior 
near these points, we follow the discussion by Fu [26] . 
For large charging energy Eq ^ F^, we can restrict the 
analysis to two charge states, for example near Ng = 1 / 2 , 
where Vp ± Ap is either 0 or 1. The charge on the dot 
can then be represented by a fermionic degree of freedom: 
A = Al ± Ar, = /!/, such that 

(AL±AR-Ag)2 = /t/(l-2Ag)-A2. (32) 

We can also express the tunneling Hamiltonian in terms 
of the fermion /. This is done by observing that 
[cj) ^ 7 r, A] = cj, ^ 7 r, such that the same dynamics is ob¬ 
tained if we replace c|, ^, 7 ^ by c\^f and treat c|, ^ and 

/ as independent fermions, because [c\ /V] = ^f. 

The Hamiltonian that governs the system projected onto 
the two charge states, 0 and 1 , can thus be written as 

Hp= Y [Hei,r + Hp^TA+Ecff{l-2Ng), (33) 

r^L,R 


with 


Ap,T.r =trY ± cIj'^ . (34) 

k 

Likewise, the current operators becomes 

It = -ietr Y - c\ j^ . (35) 

k 

The projected Hamiltonian describes a resonant level, /, 
coupled to electron reservoirs, r = L and r = R. The 
conductance is then simply given by the conductance of 
a resonant level 

r = ^FlFr 

^ ~ h AEl{l-2NgY + 2 TiXA ^ 

where F^ = F^ ± F^. Eq. ( [M| ) agrees with the pertur¬ 
bative cotunneling result in Eq. © if only cotunneling 
via the A = 1 (and not A = — 1) charge is included, i.e. 
keeping only the Eq term in ( |31| ). Parameterizing the 
L/R asymmetry as sin0p = Fp/F and cos0p = Fr/F, we 
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can write the conductance in the charge-projected basis 
as 


Gp = 


2 sin 29p 


h 1 -h sin 29p + {2Ec{l - 2Ng)/Vy ' 


(37) 


The non-perturbative result in (36) was derived un¬ 


der the assumption that the charge fluctuations are sup¬ 
pressed to the two charge states closest in energy, or in 
other words, Ec ^ r^. In order to investigate the op¬ 
posite limit where contacts have large transparency, we 
describe the Majorana junction in terms of a dual rep¬ 
resentation, which is an expansion around the perfect 
Andreev fixed point [12j . In the dual representation, the 
action of the two Majorana junctions reads [T^ 


solution by Furusaki and Matveev |53], but while they 
mapped two QPC coupled to a metallic dot to a Ma¬ 
jorana junction, we here map two Majorana junctions 
connected in series via a metallic dot to a single QPC. 
The backscattering matrix element in the QPC model 
translates to I4s = 27raAbs- Moreover, the current oper¬ 
ator is / = —idT^-12'K for both the original model and 
for the QPC model. Therefore, we can directly obtain 
the conductance for the dual model from the QPC result 
as 


Gdual — 


l + |'Fb,/2uF|2’ 


(43) 


and hence 


‘S'duaJ^L, 4 >r] = 


27tI3 


\^n I l^r (^^n) | 


ILJn 

r=L,R 


+ X! / dr cos2$r(T)-h Sc, 

(38) 


r=L,R 


with 


>S'c[‘&L,‘I’ r] “ y dr -I- 


4>T: 


-iV„ 


(39) 


where the fields 4)^ give the charge that passes through 
junction r. It is natural to introduce a total-charge field 
$+ = and a difference field $_ = $l, which 

is related to current via I = —idr^-l2TT. For Ec ^ A^. 
and at low energies |a;„| <C Eq the mode $+ gets pinned 
at irNg. We can therefore integrate out the total-charge 
field by replacing the cosine terms in the action (381 by 


their averages over $+. This is valid, when Gc, A ^ 
where A is the high-energy cut-off parameter. This leaves 
the following effective low-energy action for the difference 
field 

iuJn 

Ty^'dr (Ab,e*^-M+c.c.) (40) 


where 


^bs — + ^R, 




(«) 


The expectation value (e*^+)p appears after integration 
over $+ and it becomes 


(e*''+)o = exp(-^^ 


1 


7r|w„| -I- AEc 


AEq 
ttA ’ 


(42) 


where A = wp/a is the cut-off energy. The model (40) is 
now solvable because it maps to a single non-interacting 
QPC [23]. The spirit of this mapping is equivalent to the 


Gdual — 


/I l + (8Ac/A2)2(A2^ + 2AL.b.AH,ba 


cos 2ttNs 


(44) 


where A^^ = A^ bs + A^bs- AVe also define an de¬ 
vice asymmetry angle for the weak backscattering limit: 
sin^duai ~ Aib,bs/Abs and cos^duai ~ Ab^bs/Abs; so that 


Gdual — T 


1 


h l + (8F;cAba/A2)2(l. 


■ sin 20duai cos 2TTNg) ’ 
(45) 


We see that in both the weak tunneling and in the weak 
backscattering formulations the conductance peaks at the 
charge neutrality points and reaches e^//i for a symmet¬ 
rically coupled dot, agreeing with the expectation for the 
series resistance of two resistors each having resistance 
/i/2e^. However, in the general case they have different 
lineshapes. Fig. shows the resulting conductance for 
the charge projected model and the dual model. 

It is interesting to note that the conductance for the 
metallic dot coupled to two MBSs does not depend on 
temperature, even though it is an inelastic process. This 
is in contrast to the usual inelastic cotunneling in metallic 
dots, where the conductance is proportional to tempera¬ 
ture squared, even in the strong-coupling limit |29j . 


V. SUMMARY AND CONCLUSIONS 

Tunneling spectroscopy is an important diagnostic tool 
for determining the presence and properties of MBSs. In 
this paper, we have studied the influence of interactions 
with the surrounding circuit on the measured current- 
voltage characteristics. Confirming the results of Ref. [201 
we obtain a transition from perfect Andreev reflection 
to an isolating state of the normal-metal-MBS junction 
which can be controlled by the resistance of the circuit. 
When the environment resistance exceeds /i/2e^, the lin¬ 
ear conductance goes to zero at low temperatures. We 
studied two types of environments, namely one that can 
be described by a linear impedance and a quantum point 
contact, which allows for direct tuning of the transition. 
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FIG. 2. The solutions for the conductance of a metallic island 
connected in the two MBS. The top curve is the result for the 
weak backscattering formulation (Eq. ( |45[ )) for a symmetric 
device with i?cAbs/A^ = 0.1, while the lower curve is the 
result for weak tunneling formulation (Eq. ( |37[ )) with T /Ec = 
0.1, also for a symmetric device. For the weak tunneling case 
we have added two curves, one centered at = 1/2 and one 
centered Ng = 3/2, corresponding to projection to = 0,1 
and N = 1,2, respectively. 

Moreover, we introduced an interesting setup where 
two MBSs are coupled to a metallic island. This sys¬ 
tem has the property that the linear cotunneling current 


is temperature independent, which is very different from 
the usual cotunneling current through a metallic dot cou¬ 
pled to two metallic electrodes. In the latter case, the co¬ 
tunneling current is proportional to temperature squared 
due to the phase space of final states with an electron- 
hole pair created at each contact. In contrast, with two 
MBS junctions connecting the dot, the integral over final 
states involving a single electron-hole pair also involves 
the divergent resonant Andreev reflections and, as a re¬ 
sult, produces a temperature independent conductance. 
The metallic dot case is related to Coulomb blockade of a 
topological superconducting island considered by Fu [22] 
and it can be solved in a similar way after projection 
onto two charge states. Furthermore, the temperature 
independent result turns out to be linked to a mapping 
to an effective low-energy model consisting of a ID non¬ 
interacting fermion model with a single impurity. This 
mapping allowed us to derive a non-perturbative result 
for the low-temperature Coulomb blockade traces as a 
function of gate voltage. Finally, it is interesting to note 
that if the metallic island in this setup has more con¬ 
necting MBSs , the system is an ’’inverse out” version of 
the topological Kondo effect system recently analyzed by 
several authors 
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